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Abstract 

We give a new proof for the area law for general ID gapped systems, which exponentially 
improves Hastings' famous result [1]. Specifically, we show that for a chain of d-dimensional 
CO spins, governed by a ID local Ifamiltonian with a spectral gap e > 0, the entanglement en- 

tropy of the ground state with respect to any cut in the chain is upper bounded by '' ). 
Our approach uses the framework of Refs. |2l |3] to construct a Chebyshev-based AGSP (Ap- 
^ proximate Ground Space Projection) with favorable factors. However, our construction uses 

d the Hamiltonian directly, instead of using the Detectability lemma, which allows us to work 

' ' with general (frustrated) Hamiltonians, as well as slightly improving the 1/e dependence of the 

bound in Ref. |^. To achieve that, we establish a new, "random- walk like", bound on the 
I I entanglement rank of an arbitrary power of a ID Hamiltonian, which might be of independent 

interest: ER{H^) < (^d)*^*^'. Finally, treating d as a constant, our AGSP shows that the 
' I ' ground state is well approximated by a matrix product state with a sublinear bond dimension 

^ B = e'^^'°^ \ Using this in conjunction with known dynamical programing algorithms, 

yields an algorithm for a l/poly(n) approximation of the ground energy with a subexponential 
running time T < exp (eOd^s^^* "/^^^")) . 



'T"' 1 Introduction 

Understanding the structure and complexity of ground states of local Hamiltonians is one of the 
central problems in Condensed Matter Physics and Quantum Complexity Theory. In gapped sys- 
tems, a remarkably general conjecture about the structure of ground states, The Area Law, bounds 
^ the entanglement that such states can exhibit. Specifically, for any subset S of particles, it bounds 

CO the entanglement entropy of ps, the reduced density matrix of the ground state restricted to S, by 

" . I the surface area of 5, i.e., the number of local interactions between S and S @]. 

^ ^ Although the general area law remains an open conjecture, a lot of progress has been made on 

proving it for ID systems. The breakthrough came with Hastings' result Jj, which shows that the 
^ entanglement entropy across a cut for a ID system is a constant independent of n, the number of 

particles in the system, and scales as e'^^~J~\ where d is the dimension of each particle and e is 
the spectral gap. This result implies that the ground state of a gapped ID Hamiltonian can be 
approximated in the complexity class NP. 
In this paper, we: 
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• Give an exponential improvement to 0{ ^°^^ in the bound of entanglement entropy for the 
general (frustrated) ID Hamiltonians. The dependence on the gap even improves the previous 
best bound for frustration free ID Hamiltonians and may possibly be tight to within log factors. 

• Prove the existence of sublinear bond dimension Matrix Product State approximations of 
ground states for general ID Hamiltonians. This implies a subexponential time algorithm for 
finding such states thus providing evidence that this task is not NP-hard. 

We also establish the following properties of local Hamiltonians which may be of independent interest: 

• "Random walk like" behavior of entanglement: for a ID Hamiltonian H, the Entanglement 
Rank (ER) of is bounded by (^d)O(^). 

• Let H' be the Hamiltonian consisting only of terms acting on a subset S of particles. Then the 
ground state of H has an exponentially small amount of norm in the "high" energy spectrum 
of H': the total norm with energy above t is 2-"(*-l^^l) where \dS\ is the size of the boundary 
between S and S. 

The work here has its origins in the combinatorial approach of ^ , which used the Detectability 
lemma, introduced earlier in [S], to give a very different proof of Hastings' result for the special 
case of frustration-free Hamiltonians. The results there were greatly strengthened in [2; and [3], 
which introduced Chebyshev polynomials in conjunction with the detectability lemma to construct 
very strong AGSPs (approximate ground state projectors), leading to an exponential improvement 
of Hastings' bound in the frustration-free case to 0{{^-^^)^). 

The starting point for our results is to consider a more general situation where the Hamiltonian 
obeys the ID constraint only in a small neighborhood of s particles around the cut in question (see 
Fig.[T]). The particles to the left and to the right of this small neighborhood are acted upon by multi- 
particle Hamiltonians and Hj^ respectively. Constructing an AGSP for the new Hamiltonian is 
now much simpler, since the Hamiltonian has small norm: the AGSP is just a suitable Chebyshev 
polynomial of the Hamiltonian. In the frustration-free case, the new Hamiltonian has the same 
ground state as the original Hamiltonian, and this leads to a much simpler (and slightly stronger) 
proof of the Area Law. In the general, frustrated case, there is a tradeoff between the norm of the 
new Hamiltonian and how close its ground state is to that of the original Hamiltonian. To establish 
an area law, we must now consider a sequence of Hamiltonians whose ground states converge to the 
ground state of the original Hamiltonian, and derive an entropy bound from the tradeoff between 
the rate of convergence and the rate of increase of entanglement rank. 
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s particles 




Figure 1: The ID setting. We focus on a segment of s particles around the cut, denoting the 
muhiparticle Hamihonians to the left and right of the segment by and Hji respectively. 

2 Background: Approximate Ground State Projectors and 
their consequences 

The overall strategy is to start with a product state |?/;) and repeatedly apply some operator K 
such that p^=jp^ if-' I V') approximates the ground state and the entanglement rank of K^\ijj) is not 
too large. This property of an operator K is captured in the following definition of an approximate 
grounds state projection (AGSP): 

Definition 2.1 (An Approximate Ground-Space Projection (AGSP)) 

Consider a local Hamiltonian system H — Hi on a ID chain, together with a cut between particles 
i* and i* + 1 that bi-partitions the system. We say that an operator K is a (D, A) -Approximate 
Ground Space Projection (with respect to the cut) if the following holds: 

• Ground space invariance: for any ground state \r), K\T) = \T). 

• Shrinking: for any state \r^) E , also K\T^) e , and \\K\T^)\f < A. 

• Entanglement: the entanglement rank of K , as an element of the tensor product of two 
operator spaces (for the first and the second part of the system), is at most D. 

(The last condition implies that the operator K changes the entanglement rank of an arbitrary 
quantum state \(j)) at most by factor of D, i.e. ER(_ft'|0)) < D ■ ER((/)).j 

The parameters A and D capture the tradeoff between the rate of movement towards the ground 
state and the amount of entanglement that applying the operator K incurs. In it was shown 

that a favorable tradeoff gives an area law: 

Lemma 2.2 // there exists an {D, A)-AGSP with D ■ A < i, then there is a product state whose 
overlap with the ground state is /i — |(r|0)| > l/\/2D. 

Lemma 2.3 If there exists a product state whose overlap with the ground state is at least fi, together 
with a {D, A)-AGSP, then the entanglement entropy of \T) is bounded by 



Combined, the above two lemmas give conditions for an area law: 

Corollary 2.4 (Area Law) // there exists an {D, A)-AGSP such that D ■ A < ^, the ground state 
entropy is bounded by: 

S<0{l)-\ogD. (2) 

3 Overview 

The results here rely on the construction of a suitable AGSP that allows the application of Corollary 
|2.4[ The first critical step is to exchange local structure far from the cut for a valuable reduction 
in the norm of the Hamiltonian. To do this, we isolate a neighborhood of s + 1 particles around 
the cut in question, and then separately truncate the sum of the terms to the left and to the right 
of these s + 1 particles. Specifically, we define the truncation of an operator as follows: 

Definition 3.1 (Truncation) For any self-adjoint operator A, form A-*, the truncation of A, by 
keeping the eigenvectors the same, keeping the eigenvalues below < t the same, and replacing any 
eigenvalue > t with t. 

We then define = iT,t<i Hi)-^ + Hi -\ \- Hs + iT,i>s Ht)-\ where the s middle terms act 

on the the isolated string of s + 1 particles around the cut. The result is a Hamiltonian H that is 
now norm bounded by w = s + 2< acting on n particles with the following structure: 

H ^ =Hl + Hi+H2 + ---+H, + Hr, (3) 

where each Hi are norm bounded by 1 and acts locally on particles m + i and m + i + 1, Hl acts on 
particles 1, . . . , m and Hr acts on particles m + s + 1, . . . , n. We are interested in the entanglement 
entropy across the cut in the middle, i.e., between particles m + s/2 and m + s/2 + 1. In the 
frustration free case, it is clear that the ground state of if'*' is the same as that of the original 
Hamiltonian and it can be shown that the spectral gap is preserved for some constant value of t. 
For the frustrated case, the ground state of if'*-' is no longer that of the original Hamiltonian and 
a limiting argument (see below) will be needed to complete the proof. 

Having reduced the problem to a Hamiltonian with bounded norm u of the form ([s]) , we turn to 
the next critical step of constructing the AGSP, the use of Chebyshev polynomials to approxi- 
mate the projection onto the ground state. We begin with a suitably modified Chebyshev polynomial 
Ci{x) of degree £ with the properties that Ci{0) = 1 and |C^(x)| < e ^(^V^/") for e < x < u. The 
AGSP is then K = Ci{H) and it is clear that A = e""(^^^) . 

Bounding the ER for K requires important new ideas. We may take the approach of [H [3] as a 
starting point and expand H^ into terms of the form Hj-^ ■ ■ ■ Hj^ . For each such term, there is some 
i such that Hi occurs at most £/s times. Thus, the entanglement rank of the given term across cut 
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i is less or equal to d^^/*. The ER across the middle cut is at most times times greater, which 
gives an upper bound cP^^'^^'^. The difficulty is that the number of terms, (s + 2)^, is too large. To 
address this issue, we introduce formal commuting variables Zi and consider the polynomial 

P{Z) - {HlZo + H,Z, + ■■■+ HrZ, + ,Y = /ao,...,a. + i^o"^r • ■ ■ ^s+t' ' 

aoH i-as + i=l 

In particular, — J2ao a^+i /oo,--.,as+i ■ This expression has fewer terms, namely, (^^^^^)- As 
before, for each multi-index (oq, . . . , cis+i) there is some i such that flj < l/s. If i is fixed, a linear 
combination of the corresponding operators fag.. ..,(1^+1 can be generated as follows. We restrict our 
attention to only those terms in P{Z) where Zi appears at most £/s times and assign arbitrary values 
to the variables Zq, . . . , Zs+i- The ER of the resulting operator is estimated using the representation 
P{Z) = {A + HiZi + BY, where A and B commute. We then use a polynomial interpolation 
argument to express each fag,. ..,0^+1, their sum H^, and finally, the operator K. Thus we prove that 
the ER of K is at most D = 



Applying Theorem 2.4 to the above AGSP with £ = 0{s'^), s = 0(log^(d)/e) yields our Area Law 
for frustration free Hamiltonians, providing an entanglement entropy bound of the form 0(log'^ {d) /e) . 

To address the frustrated case, a third critical result is needed: that the ground states of H^-*^ are 
very good approximations of the ground state of the original Hamiltonian. Intuitively, the structure 
of the small eigenvectors and eigenvalues of i/^*) should approach those of H as t grows and we 
show that to be the case, showing a robustness theorem; that the ground states of ff^*-' and H 
are exponentially close in t and the spectral gaps are of the same order. 



We would like to apply Theorem 2.4 to an AGSP for H^^\ for t sufficiently large, however, if we 
try to do this in one step, the ER cost becomes a large function of t. Instead we use a well chosen 
arithmetic sequence to,ti, . . . and the associated AGSP's to iJ^*') to guide the movement towards 
the ground state. The robustness theorem allows for very rapid convergence, the result of which is 
the area law in the general (frustrated) case. 



4 Approximate Ground State Projector 

Consider a Hamiltonian H acting on n particles with the following structure: H = Hl + Hi + H2 + 
■ ■ ■ + Hs + -fffl,, where Hi acts locally on particles m + i and m + i + 1, Hl acts on particles 1, . . . , m 
and Hji acts on particles m + s-l- 1, . . . , n. Assume that H has a unique ground state |r) with energy 
£9 and that the other eigenvalues belong to the interval [ei , u] . Let e = ei — Cq denote the spectral 
gap. We wish to bound the entanglement entropy of \T) across the middle cut, i = s/2. (In our 
notation, cut i separates the particles m + i and m + i + 1.) 

We define the AGSP a.s K ^ Ce{H), where Ci is a polynomial that satisfies the conditions below 
for a suitable value of A. 

1. Q(eo) = 1; 
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2. \Ct{x)\ < \/A for ei < a; < u. 

It follows that K\T) = \T) and that the restriction of K to the orthogonal complement of |r) has 
norm less or equal to \/A. 



Lemma 4.1 There exists a degree i polynomial Ce that satisfies the above conditions for 

= 2 e~^^V (<^i-<^o)/ ("-eo) _ 

Proof: We construct by a linear rcscaling of the Chebyshev polynomial Ti, which is defined 
by the equation Ti(cos9) = coii{£d). It follows immediately that |T^(.t)| < 1 for a; € [—1,1]. If 
a: > 1, the equation cos 6* = x has a complex solution, 9 = it, where cosht = x. In this case, 
Ti{x) = cosh(£t) > ie". Since t > 2tanh(f/2) = 2^(x - l)/(a; + 1), we conclude that 



Now, let 



Te{x) > ^e2V(^-i)/(^+i). 



Ciiy) = XV ^l^ere / y = — . 



The function / maps ei to 1 and u to —1, hence |Q(y)| < T^(/(eo)) ^'^^ ^ ^ [^Ij^]- The bound for 
Ti{x) with X = /(eo) matches the expression for A because — -^^^i^^-^^^"^ — 



rK/(eo)) 

si-ep _ ^ _ 

u-eo /(m)-/(<:o) x+1 ' 



Lemma 4.2 T/ie entanglement rank of K = C^{H) (where Ci is an arbitrary degree £ polynomial) 
is bounded by D = (ci£)0(max{^/s,v^})_ 

Proof: W.l.o.g. we may assume that s < \/£. If that is not the case, we can reduce s to Vi by 
joining some of the Hj^s with either or Hn. This does not change the actual entanglement rank 
or the required bound. After this reduction, the bound can be written as {dl)^'^^/^\ 

K = Ce{H) is a linear combination of i + 1 powers of H, and we will bound the entanglement 
rank added by each, focusing on the worst case H^. Let us first consider the expansion = 
je ' ■ ' ^^i' '^^'^y tsrms to be useful, but we can group them by the number of 

occurrences of each Hj. To this end, we introduce a generating function, which is a polynomial in 
formal commuting variables Zq, . . . , Zg+i: 

Pe{Z) = {HlZo + H,Z, + ■■■ + HnZ,+,f = ^ /ao,...,a.+.^o°^r • • • 

aoH \-as+i=( 

Each coefficient /ao,...,os+i is the sum of products Hj^ . . . Hj^, where each Hj occurs exactly Oj times. 
We are interested in estimating the ER of = ^^^^ fao,-.,as+i- 

We start by noticing that for each multi-index (ao, . . . , ag+i), there is some i G {1, . . . , s} such 
that Oj < l/s. Thus, = J^^^i J2k=o Qi,ek, where Qi^ek includes some of the operators fao,...,as+i 
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such that ai — k and X^j^ti % ~ ^~ ^- (This decomposition of is not unique.) We will, eventually, 
bound the ER of each Qi/k- To do that, we first define a generating function that includes all the 
matching /ao,...,a,+i 's: 

Pi,ikiZ) ^ ^ /ao....,a, + i JJ-^"'- 

ai—k j^i 
T,j^i aj=i-k 

This sum has t = (^"s^'*) terms. The variable Zi is excluded, or we may consider it equal to 1. 
When the remaining variables are assigned definite values, Z e C*^^, we obtain a linear combination 
of the operators /ao,...,aa+i ■ The key observation is that such linear combinations have full rank, i.e. 
there are t distinct values of Z e such that the corresponding {Pi^(k(Z)} form a basis in the 

space of operators of the form ^ao,.. .,0^+1/00,... ,0^+1 7 where Cao,....as+i G C and the sum runs over 
the support of Pijk- In particular, Qijk is a linear combination of t operators of the form Pijk{Z). 

For a fixed Z, the operator Pi ik{Z) can be obtained as follows. We write Pe{Z) = {A + Hi + BY , 
where A = X]j<i ^j^j and B — X]j>i ^j^j, and then collect the terms with Hi appearing exactly 
k times. Since A and B commute, such terms have the form A"'°B''''Hi ■ ■ ■ HiA"''' B^'' . There are 
(2^fe'+i) distinct terms like that, and the ER of each term across cut i is at most cP^ . The ER across 
the middle cut is bounded by that number times (d^)''^"/^' < . Combining all factors, we find 
that 

mQ^M) <{^']^'^ {2k + ^ - < {Mf^'/^\ 

Here we have used the fact that k < t/s and s < y/I. The summation over i and k does not change 
this asymptotic form. ■ 



Lemma [4^2] gives a non-trivial tradeoff between the entanglement rank D and shrinking coefficient 
A of the operator K. By suitable choice of parameters this will give the desired {D, A)-AGSP such 



that D ■ A < ^ and Corollary 2.4 will apply. One issue that we will have to address is the bound 
t on the norms of Hl and Hij. We first tackle the case of frustration free Hamiltonians, where we 
can assume W.L.O.G. that t = 0(1) = \Hl\ = \Hr\: 

Let H' — ^ Hi he a, frustration free Hamiltonian with spectral gap e. For t chosen in a moment, 
define Hl — {J2i<m ^i)^'^ Hjj, = {'l2i>m+s+i ^i)^*' to be the truncation of the Hamiltonian 
acting on the left and right ends of the line. Set H = Hl + Hi + H2 + ■ ■ ■ + Hg + Hj^ so it is in the 
form as above. Clearly H has the same ground state as H' . Since eo — 0, Lemma [GJ] (below) yields 
that for t = 0(i) the Hamiltonian H has a gap that is at least a constant times e. 

Theorem 4.3 For a frustration free Hamiltonian H' — ^Hi with gap e the entanglement entropy 
IS O(^). 



Proof: From Lemma 



6.1 



for t = O(-), we have that H has the same ground state and gap of the 



same order as H' . Recall lemma 



4.2 



applied to H describe an AGSP with bounds A = e vs+2t and 
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fl/2.,3/2 



and 



D = {I + + Set I = sV2, c = y^/^/^T2^ so that A = 

Write the condition DA < 1/2 as logD < log - 1. log ^ = ce^/^s^/^, and logD = 0{s{logd + 
log(s2 + s))). Thus we can satisfy the condition with s ^ O^-^^), and therefore logZ? = 0{C-^^)) 



and the result follows directly from Corollary 2.4 



5 Low bond dimension MPS for frustration free ID Hamil- 
tonians 

We can use these results to show the existence of a matrix product state of sub-linear bond dimension 

1 2 

of size exp(0(e~3 log^ n)), that approximates a ground state |r) of a gapped frustration free ID 
Hamiltonian to within pgiyj^^) ■ show the existence of a matrix product state of bond dimension 
B within 5 of the |r), it suffices to show the existence of a state of entanglement rank B within - 

of|r). 

We've shown the existence of a state with constant overlap with |r) and entanglement rank 
0(]2E_A), To this state, we would like to apply an AGSP with A = ^i^i^) ■ J^^t as in the proof of 



4.3 



we choose I = and Lemma 



4.2 



poly(n) ■ 

estabhshes that a A = e"'^^ ^ , D = exp{s log d) 



Theorem 

AGSP exists. Setting s = 0{e^^ logs n), we have A — poiy(n) ^i^h D — exp(0(e 3 log^ n)) 



6 Frustrated Case 
6.1 The operator H^^^ 

We consider the ground state |r) of a local Hamiltonian H' = '}2,H[, where H[ acts locally on the 
particles i and i + 1, and < H[ < 1. We assume the Hamiltonian H' has a unique ground state 
with energy eg and next lowest energy ei; let e = ei — eo denote the spectral gap. It is easy to see 
that in such a case e < 1. 

Ideally we wish to replace H' with some Hamiltonian H with the same ground state |r) and 
spectral gap e, but with smaller norm, so that the AGSP from Section |4] yields a good bound on 
the entanglement entropy of |r). Towards that goal we consider Hamiltonians of the following more 
general form: 

H = Hl + H,+H2 + --- + H,+Hr, (4) 

where Hi acts locally on particles m + i and m + i + l, acts on particles 1, . . . , to and Hji acts on 
particles to + s + 1, . . . , n. We further require that Hl and are positive and \Hi + . . . iJs-2| < s. 



8 



We now show that |r) is the ground state of such a Hamiltonian H with the added properties 
that the ground energy of Hl, Hr and YlilZt Hi are all and < Hi < 1 for i = 1, 2, 3, s — 2, s — 1, s. 
By enforcing these properties, we have eo < 6. We do this by setting: 

• Hl = J2iLi H'i — cl where c is the ground energy of X]"=i ^l, 

• Hr =^ Er=m+s+i H[ - c'l where c' is the ground energy of Y.l=rn+s+i ^1 

• Hi = Hl^_^i for the six values i — 1, 2, 3, s — 2, s — 1, s, 

• Hi — iJ„i+i — for 4 < i < s — 3, where d is the ground energy of X]i=4 ^i- 

It is easily verified that H is of the form Q, and since the difference between H and H' is a 
multiple of the identity, H has ground state |r). The Hamiltonian H has ground energy < 6 since 
the tensor product of the ground states for the disjoint operators Hl, J^tZl ^i, and Hji only can 
have non-zero energy on Hi + H2 + H3 + Hs-2 + -ff,s-i + Hs- 

To bound the norm (so as to effectively apply the AGSP from Section |4]) we use the previously 



defined notion of truncation (Definition 3.11. If for any self-adjoint operator we let Pt be the 



projection into the subspace of eigenvectors of A with eigenvalues < t, then 

A^' ^PtAPt + t{\-Pt) . (5) 

Define = {Hl + + H2 + ■ ■ ■ + Hs-i + {Hs + i?^)^*; it is easy to verify for t > 0, 

jjit) < fj and iJ*^*^ < {2t + s)l. Unfortunately, the truncated Hamiltonian i/*^*' no longer has the 
same ground state |r). Intuitively, the structure of the small eigenvectors and eigenvalues of 7?'*' 



should approach those of if as i grows. The Robustness Theorem (Theorem 6.1), stated below and 



proved in Section 6.3 verifies this intuition, showing that for t bigger than some constant, the gap 
of is of the same order as the gap of H and the ground states of i/'-*^ and H are exponentially 
close in t: 

Theorem 6.1 (Robustness Theorem) Let |r),eo,ei be the ground state, the ground energy and 
the first excited level of H, and let |0),eg,e']^ be the equivalent quantities of H^^\ Then for t > 
0(^_(^o_ _|. 1)) have 

a. e'l - ef, > 0(ei - eq) 

b. |||0)-|r)f <2-<^(*). 

6.2 General ID Area Law 



We would like to apply lemma 4.2 to for t sufficiently large, however, if we try to do this in one 
step, the entanglement rank cost becomes a large function of t. The overall plan is therefore to use a 
well chosen sequence tQ,ti, . . . to guide the movement towards the ground state. More concretely we 
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define a sequence of states |V'o)i IV'i)! ■ • • that converge to |r), while carefully controlling the tradeoff 
between increase in entanglement rank and increase in overlap with |r). 

Denote by 1 — fii, the overlap between the ground state |(/)*') of and |r). We use the AGSP 



K — Ci{H^*^^^) from lemma 4.2 to move from state to jV'i), where lipi) has overlap at least 

1 — i-ii with 10*'). We will show that the increase in entanglement rank of each move is small enough 
to bound the entanglement entropy of the limiting state which is the ground state of H . 
We now put all the ingredients together to prove an area law for general ID systems: 

Theorem 6.2 For any Hamiltonian of the form H = + Hi + • • • Hg + Hj^ with a spectral gap of 

1 ^ d 

e, the entanglement entropy of the ground state across the (s/2,s/2 + l) cutis bounded byO{^^^). 
We begin with a lemma: 

Lemma 6.3 There are constants tg and c and states \ipi) with entanglement rank Ri, i = 0, 1, 2 . . . 
satisfying: 

1- > 1-0(2-*), z>0. 

2. logi?o = OC^), logi?, = logi?o + 0(E-=i^,logd), withtj^O{^[t,+s)/e). 
Proof of Lemma: 

We begin by choosing constants to, c such that Q.{tQ + ic) > z + 4, for the fl{t) appearing in 
Theorem |6.1| Setting ti = to + ic, we therefore have 

|||r)-|0*of <2-(^+*) (6) 

for all i. Similar to the frustration free case, since to is constant, choosing i = s^ and s = 0( '°^^ '^ ) 



gives DA < i in Lemma 4.2 and thus by Lemma 2.2 there exists a product state |V') such that 

this time with I chosen so that 



4.2 



KrilV')! < ^ where logD = 0((^^)). Returning to Lemma 

A = e V=+2*o-=o ^ we estabhsh that the state IV'o) = (H<!°)|vi)|| property that 

1110*") - I0'o)||^ < while having entanglement rank i?o with logi?o = OQ^^). 

We now inductively define |'(/'i) from and show that 1110**) — |V'i)ll^ ^ 2^^^'^. Applying the 

triangle inequality to the induction hypothesis |||0*'~O ~ iV'i-i)!!^ ^ 2~^^^ along with the already 
established proximity of |0*'-i), |0**) to jP) of yields 

1110*0 -|^.-i)ll'<2-'-^ (7) 
Our goal is, with only a small amount of added entanglement, to move jf/'i-i) a little bit cl oser 



to 10*') which we will accomplish by using a well chosen AGSP. With ti = 0(y ^^i^), Lemma 4.2 
establishes the existence of a (I?, ^) AGSP K for iJ^*') with a loose bound of logZ? < 0(£ilog(i); 
we apply this AGSP K to move from |'0i-i) to |?/'i) by setting |-0i) = . The shrinking 

property of the AGSP along with (m) establishes || 10*') - IV'i) 11^ < 2 
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All told we have generated states Itpi) with entanglement rank i?,; = _Ro + J2j<i^i^'^S^ ^^'^ 
1110*') - \tp'i)f < 2^'"^. Finally, \{iJr\r)\ > 1 - |l|r) - \^p^)f/2 > 1 - 2-' where the last inequality 
again used (|6|. 

Proof of Theorem: 

The above lemma gives a series of states of bounded entanglement rank that converge to the 



ground state |r). Thus if {A;} are the Schmidt coefficients of F, Lemma 6.3 



Y.xj>m\r)f>i-o{2-^). 

i=l 

The entropy of |F) is then upper bounded by summing 0(2"*) logi?i (i.e. the maximal entropy 
contribution of mass 0(2"') spread over logi?i terms). We arrive at a bound of the entanglement 
entropy of |F) given by: 



J2 0(2-^) (log Ro + 0(^ i, log d)) = log Ro 0(2-') + ^ 0(2"') ^ ^ ^o + ^c + s ^ 

i j = l i i j = l 



= Oilo,R,) + Oi^) = Oi'^). 
6.3 Proof of Theorem 16.11 

Before giving the proof, we use a simple Markov bound to show that in a gapped situation, a state 
with low enough energy must be close to the ground state: 

Lemma 6.4 (Markov) Let B be a self-adjoint operator with lowest two eigenvalues eg < ei; denote 
its lowest eigenvector by Given a vector \v) with low energy, i.e. such that {v\H\v) < eo + S, 
then \v) is close to lip) in the following sense: 

\m~\v)f<^^. 

ei - £0 

Proof: Write \v) — a\ip) + Vl — a'^\ip-^) where jV"^) is orthogonal to |-0). The energy of \v) then 
satisfies 

a^eo + (1 - a^)ei < {v\H\v) < eo + S, 

and thus (1 — a^) < j^^r^- The result follows from noting that — |w)||^ = (1 ~ a)^ + (1 — o,^) = 
2 - 2a < 2(1 - a^). ■ 

Proof of Theorem I6.lt Define A to be the sum of the two terms A = H2 + -ff^-i- Notice that 
the operators: 

{(fr(*) - A),H^ + Hi,Hs + H^, H-A:t>0} 
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all commute with each other and therefore all the operators are simultaneously diagonalizable, i.e. 
they have a common collection of eigenstates. We fix Pt to be the projection onto the subspace 
spanned by those eigenstates oi H — A with eigenvalues less than t. This collection of eigenstates 
clearly have eigenvalues less than t for the operators + Hi and Hg + and therefore + Hi = 
{H^ + Hi)-* and Hg + H^ = {Hs+ H^)-* on the range of Pt. This allows the important observation 
that 

jjWp^ = HPt. (8) 

Proving part [a] 

The main idea is to consider the normalized projection of the two lowest eigenvectors of i?^*^ 
onto the lower part of the spectrum of H — A using Pt. Under the assumption that t is sufficiently 
large and the gap of iJ is sufficiently small (relative to the gap of H) , we show contradictory facts 
about these normalized projections: that they are simultaneously far apart from each other (because 
applying Pt did not move either very much) and close to |r) (because they both have energy with 
respect to H that is close to eg). This contradiction allows us to conclude that for t sufficiently 
large, the gap of -ff must be of the order of the gap of H. 

For every normalized state \v), we define \vt) = Pt\v) and \vh) — (1 — Pt)\v) . Note that by ([s]), 

{vt\H\vt) = {vt\H^'^\vt) (9) 

Our main technical tool is the following lemma that connects the energy of a state \v) to that of 
\vt): 

Lemma 6.5 For any state \v) , 

1- \\W)\\<\f^^, 

2. {vt\H\vt) < {v\H\v)+0{^/^^^^p^). 
Proof: The first result follows from 

(Hi/(*V> > {v\iH^'^ - A)\v) = (i;,|(i/W - A)\vt) + {v,^{H^'^ - > («,|(i/(*) - A)W) > t\\v,^ 

the last inequality by the definition of Pt. 
For the second result, 

{vt\H\vt) = {vt\H'^''^\vt) < {v\H^'^\v)+2\{vt\H^'^\vH% 

the first equality from Q and the second inequality from writing \vt) — \v) — \vh) and expanding. 
We now bound the second term on the right hand side. Notice that (ut|i7'^*)|w/i) = (wtji/*^*'' — 
A\vh) + {vt\A\vh) ^ {vt\A\v h). To \{vt\A\vh)\^ we apply Cauchy- Schwartz to get \{vt\A\vh)\^ < 
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To prove part [a[ assume e'^ — ef) < jq^^i — co) and denote by \4>^) the eigenvector of i?^*^ with 
eigenvalue €[. Write \(j)t) — Pt\<j)), \(t)}) ~ Pt\(f>^). Lemma 



6.5 



estabhshes 



mt)<eo + 0{J^) 
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t 



Setting \v) = 



j^'jll , = iij^i^ii , and using the above in Lemma 



6.4 



yields 



\T)-\v}f<0{l)^^^ 
(\ - eo \Jt 



llir) - \v')\f < ^ + 0(1) ^ /° 

This establishes, for sufficiently large t = ), that and are both near |r) 

contradicting the fact that they are also almost orthogonal. ■ 

Proving part [b} 

We are interested in showing that the ground states of H and iJ^*-' are very close together. 
Clearly, the ground states of the nearby Hamiltonians H — A and 77^*^ — A are identical since they 
only differ among the eigenvectors with values above t and so the question becomes how much the 
addition of A can change things. This reduces to how much the operator A mixes the low and high 
spectral subspaces oi H — A (i.e. how big the off-diagonal contribution of A is when it is viewed 
in a basis that diagonalizes H — A). The core component of the argument will be the Truncation 



Lemma (Lemma 6.7): that the ground state |r) is exponentially close to the range of Pt (i.e. the 
low spectral subspace of H — A). We will combine this result with the fact that H and iJ^*' are 
identical on the range of Pt to argue that ground states for H and are exponentially close. 
The following lemma captures the bounds necessary for proving the Truncation Lemma. 

Lemma 6.6 With H, Pt, \T) , eg as above we have the following: 

1- \\{i-Pt)m^ <'^^^^^N^^^^^, 

2. Fort>u, \\{l-Pt)HPu\\ = ||(l-Pt)^^'u|l < 2e-^. 
Proof: For 1., by definition, 
e = {T\H\T) = {T\PtHPt\T) + (r|(l - Pt)H{\ - Pt)\T) + {T\PtH{l - Pt)\T) + (r|(l - Pt)HPt\T). 
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This gives the bound 

e > e\\Pt\r)f + t\\ii - p,)\T)\f - 2|(r|(i - Pt)AP,\r)\, 

where the second term of the right hand side foUows from the inequahty (r|(l — Pt)H{l — Pt)\T) < 
(r|(l — Pt){H — A){1 — Pt)\r), and the replacement of H with A in the third term follows from 
il-Pt)iH-A)Pt = {l-Pt)Pt{H-A)^Q. Writing ||Pt|r>||' = 1 - ||(1 - Pt)|r>||' and rearranging 
terms yields statement 1. 

For statement 2., the first inequality follows simply from writing H = [H— A)+A and noting that 
H-A commutes with P„. We write (1 - Pt)APu = (1 - Pt)e-''^"-^h'''^"-'^^ Ae-''^"-'^h'''^"-^^ Pu, 
for an r > to be chosen later, and noting therefore that 



\\{l-Pt)APu\\ < ||(l-Pt)e""^''"'^^ll • ||e'-(^-^Ue-''(^-^)|| < e-''^*-") He-^^^--^) Ae-"(^-^) || . 
The Hadamard Lemma gives the expansion 

^r{H-A)^^-r(H-A) = A + r[H - A, A] + —[H - A, [H - A, A]] (10) 
+ [H-A, [H-A,A]]] + --- 

= Qo + rQi + '^Q2 + '^Q3 + --- , 

and we turn to bounding the norm of these operators Qi. 

If we expand H — A and A as the sum of its constituent local terms Hj , each Qi can be written 
as a sum of Ui terms, each a product of Hj^s; we now bound n.;. Notice that uq = 2 and that Qi-i 
consists of terms, each of which is a product of at most i Hj^s. For such a product, there are at most 
2i terms in H — A that do not commute with it. This implies the recursive bound Tij < 4mi_i and 
thus rii < 2 ■ 4*i!. Since each of the terms is norm bounded by 1, we have \\Qi\\ < 2 ■ 4*i!. Plugging 
this bound into (lOl we have e'"'^^""^^ Ae^'"'-^^"^-' < 2^j(4r)*; choosing r = ^ gives a bound of 4 for 
( 10 1 and establishes statement 2. ■ 

Lemma 6.7 (Truncation Lemma) For t > 17, 

||(i-PO|r)|| <2-"w. (11) 



Proof: We show a discrete version of ( 11 ); that there exist constants s — 16 + eg and d — 16 such 
that for integers n > 

\\il-Ps+nd)\T)\\<2-- (12) 
The result will then follow since s + nd < t < s + {n + l)d implies 

m - Pt)\r)\\ < ||(i-F,+„<i)|r)|| <2-" <2-(^-i) = 2-"(*). 
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To prove ( |12[ ), we will proceed by induction. Clearly the initial case of n = holds. Assume that 
11(1 - Ps+nd)W)\\ < 2"", for n < no. Define P[o] = Ps and Py] = Ps+jd - Ps+u-i)d for 1 < j < no; 
thus Ps+ngd — Sj=o ^b] ^'^'^ ^^"^ induction hypothesis implies 



Py]\r)\\<2-^+\ 



(13) 



for j < nQ. By Lemma |6.6| 



n/-, p N|T^\||2 ^ 2|(r|(l - Ps+nod)^^'s+nod|r)| 

s + nod - eo 

Our goal is to bound the numerator of the right hand side by 16 • 2~^""; (12) then follows since the 
denominator is at least 16. Write 



(r|(i-p,+„„rfMP,+„„,|r)| ^\Y,{m-Ps+nod)APy]\T)\ 

J=0 



< X] 11(1 - P.+„„<i)AP[,]||||(i - p,+„„,)|r)||||P[,]|r>|| < 5] 11(1 - Ps+nod)APy]\\2-^-"«+^ 

j=0 j=0 



where the last equation used (13) and the fact that ||(1 — Ps+n„d)\^)\\ < II (1 ^ Ps+(no-i)d)\^)\\ < 
2-(no-i)_ Applying Lemma |6.6| to bound the first term in the sum on the right hand side yields: 



no no 



|(r|(l -P,+„„rf)AP,+„„,|r>| < ^e-^^2-(««+^-3) = 2-2""(8)^2("«-^")(i-^B>b). 

With the choice of d = 16 In 2, the sum in the last term on the right hand side is a geometric series 
that is bounded by 2 which yields the desired bound of 16 • 2^^"° and completes the proof of ( 12 ). ■ 

We now use the Truncation Lemma to show that the the projected state Pt |r) is exponentially 
close to eigenvectors of both H and H'-^K 

Lemma 6.8 The state iFf) = ^^jp^Pt|r) is an approximate eigenvector of both H and H^^^ in the 
following sense: 

\\H\Tt) ~ eo|r,)|| = ||H(*)|r,) - eo|ri)|| < 2-"W. (14) 
Proof: We begin by writing H{\ - Pt)\T) = J^Zo H{Pt+^+l ~ Pi+,)|r) and thus 



\H{1 Pt)\T)\\ < J2 \\H{Pt+.+i - Pt+.)ll • ll^*+.|r)|| < Y.^t + (z + 1) + 2)2"^(*+') < T 



-0(t) 

i=0 1=0 



where the bound on the first term follows from the fact that H — A < (t + « + l)-l on the range of 
Pt+i+i and the bound on the second term from the Truncation Lemma. 
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We write eo|r) = H\T) = H{1 - Pt)\T) + HPt\T). We have bounded the first term on the right 
hand side by 2^^^*^*-' and it follows simply that 



|i/Pi|r)-eo|r) 



\H^'^Pt\T) 



mil < 2-"(*) 



where the first equality is from (|8|). Multiplying (|15|) by pryp^y ^ ^ 



< 



(15) 

a constant close to 



one yields ( 14 1 since the constant can be absorbed into 2 



We outline the remainder of the argument. The approximate eigenvalue property of ( 14 1 can be 
used to show that \Tt) is close to an eigenvector of with eigenvalue in the range [eo — 2^^^'*) , eo + 
2-f^(*)j. By combining a lower bound for the ground energy of H'^^^ with the fact (part a.) that iJ^*^ 
has a gap of reasonable size, we are able to show that the eigenvector of _ff near \Tt) is the ground 
state of The proximity of \Tt) to the ground states of both and H then establishes the 

result. 



dof 



We begin by setting 5 2^^^^*^ The approximate eigenvalue property of (14) implies that 
there is an eigenvalue of i?*-*) within 5 of eg, for if not, writing \Tt) — '^^Ci\vi) where \vj) are the 
eigenvectors of i?^*) with eigenvalues A^, \\H^'^'>\Tt) ~ tQ^t\^ — X]i(^j ^ '^o)c| > Yl,i ^1 — which 



contradicts (14). We now show that there is only one eigenvalue in this range and that it is in fact 



the ground energy for H'^^\ This will follow by lower bounding the energy of if*^*-* by eo + (5— (e'j^ — ej)). 
Decompose an arbitrary state = Pt\'4') + [1 — Pt)\4') = + \'^2) ■ Then 

im^'M) = mH'^''^ - m) + = (v-iK^?^*' - + min'^'^ - a)iv^2) + (v-iai^) > iii^2ir 



for which i||V'2|| > ^ will therefore have energy at least eo. In the remaining case of 



11^1 



> 1 



^ we bound the energy of jV') as follows: 



> {i'i\H^'^\^i) + {HH^'^\^2)-2\{^i\A\^P2)\ > (V'i|i?'*Vi) + (V'2|i?(*^|V^2)-4y^. (16) 
Since |V'i) is in the range of Pt, ^ implies (V'i|i?'*^|'0i) — (V'll^^lV'i) > ^ollV'll^- Combining this 



bound with (16) gives (tp\H lip) > eo — 0{y) and thus we've shown eg > eo — 0{y)- Since (from 
part a.) e'l — eo > 0(ei — eo), a choice of t = 0{ (-g^l°^-)2 ) ensures that the the ground energy of iJ^^*^ 
is at least eo + <5 — {e'l — Eq). 

We've established that both ground energy and the energy of \Tt) with respect to is in the 



interval [eo — 6, eo + 6]. Applying Lemma 6.4 gives 

lll^>-|rt)||<<5, 

and b. then follows from recalling that |||r) — |rt)|| < 5 as well. 
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7 Sub exponential algorithm for finding the ground energy 
of gapped ID Hamiltonians 



We now show that the ground state can be weh-approximated by a MPS with a subhnear bond 
dimension, and, consequently, a l/poly(7i) approximation of its ground energy can be found in a 
subexponential time. To simplify the discussion, we shall treat d as a constant. 

As discussed in the frustration free case, to show that |r) can be approximated to within ^^1]^^^) 
with an MPS of bond dimension B = 0(exp(log* n/ei)), it suffices to show for each cut + 1) 
the existence of a state with entanglement rank B across that cut that is within ^^^^^^^ of \T). 



Theorem 



6.1 



yields that for t ~ O(logn), — r|||| < ^^^^^^^ , and therefore we turn to finding 

a state with entanglement rank B a cross the cut + 1) that approximates l^*-'-*). Applying 
the AGSP K = Ci{H^*^) of Lemma 4.2 we have A = for £ = Cl(log n^^t^^)), and 



D = 0(exp(^/s + s)). The optimal choice of s = log-i n/eJ gives the desired B — e'^^^°^'^ "/^ 

We can now use this result to bound the complexity of actually finding the ground energy. For 
simplicity, we treat d,e as constants. Using recent dynamical programing results [?, ?], we infer 
that there exists an algorithm that runs in time T = (dBn)^^^ ^ < exp(e'^'-'°s^^ ' )) and finds a 
l/poly(n) approximation of the ground energy. Since e'^d^s'"' is smaller than any finite root 
of n, it follows that 

Corollary 7.1 Finding a l/poly(ri) approximation to the ground energy of a ID, nearest-neighbors 
Hamiltonian with a constant spectral gap is not HP -hard, unless 3-SAT can he solved in a sub- 
exponential time. 
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